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INTRODUCTION. 


^ in'ojectilc  is  subjected  during  the  first  few  milliseconds  of  free 
flight  to  a high  pressure  jet  flow  from  the  muzzle.  A.ssessment  of  the 
projectile's  behavior  in  this  environment  requires  an  accurate  descrip- 
tion of  the  jet  which  i.s  typically  provided  by  theoretical  calculations. 
Their  accuracy  depends  among  other  things  on  the  accuracy  of  the  assumed 
boundary  conditions,  i.e.,  of  the  conditions  in  the  flow  through  the  muz- 
zle of  the  weapon.  Therefore,  one  is  interested  to  assess  the  accuracy 
of  theoretically  predicted  muzzle  flows.  Such  detailed  predictions  are  rou- 
tinely provided  by  mathematical  models  of  the  interior  ballistics.  Ex- 
Derimental  determination  of  the  flow  is,  on  the  other  hand,  technically 
difficult  and  is  usually  limited  to  a measurement  of  the  projectile's 
velocity  and  to  pressure  measurements  in  the  gun  tube  near  the  muzzle. 

The  measurements  do  not  determine  the  gas  flow  completely.  Thus,  differ- 
ent computer  model.s  of  interior  ballistics  may  predict  muzzle  velocities 
and  pressures  which  arc  in  agreement  with  observations,  but  postulate 
at  the  banie  time  significantly  different  temperatures  and  densities. 
Independent  measurement  of  temperatures  and  densities  iu  muzzle  flows 
are  therefore  of  great  interest  theoretically  as  well  as  for  practical 
reasons.  ' ' = - - 


In  a recent  BRI.  Report  fl]  E.  M.  Sclimjdt,  ct  al,  have  described  a 
new  temperature  measuring  technique  which  is  designed  to  measure  the 
initial  temperature  in  the  gas  efflux.  The  technique  is  based  on  a mea- 
surement of  sound  speed  in  the  g.as  and  a subsequent  calculation  of  the 
temperature  from  the  sound  speed.  Ke  shall  give  in  Section  2 of  this 
report  a short  description  of  the  technique. 

Measurements  of  temperature  constitute  a new  data  source  (in  addition 
to  pressure  measurements)  for  the  validation  of  interior  ballistics  cal- 
culations. In  order  to  use  the  observations  for  this  purpose,  one  has 
to  know  their  accuracies.  In  this  report,  we  shall  investigate  the  in- 
trinsic accuracy  of  the  temperature  measurement  technique.  The  partic- 
ular goal  of  the  investigation  is  to  find  quantitative  limits  within 


E.  M.  Schmidt,  E.  F.  Gion  and  D.  1).  Shear,  "Acoustic  Thcrmometric  Mea- 
surements of  Propellant  Gas  Temperatures  in  Guns".  BRL  Report  No.  1919, 
August  1970.  . (AD  (IA030.SS9) 
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which  tlic  temperature  measurements  slioiild  agree  with  tcmjicratures  pre- 
dicted by  interior  ballistics  calculations  of  the  same  event. 

The  accuracies  of  the  temperature  measurements  are  affected  by 
three  sources  of  errors.  “ 

first,  one  has  a limited  observational  accur.icy  which  produces  a 
computed  sound  speed  with  .an  error  in.aTgiii.  This  error  propagates 
through  subsequent  calculations  and  producc.s  a corresponding  error  mar- 
gin for  the  temperature.  An  investigation  of  this  error  is  done  in  Sec- 
tion 3,  and  it  is  based  on  variance  propagation  formulas. 

Second  the  measurement  te.chnique  possesses  an  intrinsic  averaging 
property  due  to  the  finite  size  of  the  measurement  appavatu.= . A conse- 
quence of  this  property  is  that  the  recorded  tempeiamro  is  not  the  in- 
stantaneous initial  temperature  in  the  gas  offlu.x,  hut  an  average  tem- 
perature for  th''  first  few  milliseconds  of  the  flow.  Therefore,  short 
duration  temperature  excursions  (predicted  by  some  mathematical  models 
of  interior  ballistics)  might  not  be  detectable  by  the  apparatus.  A 
■quantitative  investigation  of  the  averaging  property  of  the  technique  is 
done  in  Sections  4 and  5.  In  Section  4 we  describe  typical  theoretical- 
ly calculated  interior  ballistics  flows,  and  in  Section  5 we  simulate 
and  analyze  the  re.sponse  of  the  measurement  apparatus  to  the  calculated 
flows . 


A third  source  of  error  is  the  approximate  nature  of  the  equation 
of  state  which  provides  the  relation  betwoon  sound  spe-c-d  and  temperature. 
An  invest ig.at ion  of  this  systematic  error  source  requires  an  analysis 
cf  the  jjliysical  properties  of  the  combustion  products,  and  it  will  not 
bo  pursued  here.  We  notice,  however,  that  for  a comparison  of  tempera- 
ture measurements  of  a real  firing  with  interior  ballistics  calculations 
this  error  source  is  not  important,  because  it  affects  measurements  and 
calculations  in  the  same  way.  The  calculations  depend,  on  the  otnor 
hand,  also  on  many  other  assumptions  about  the  properties  of  the  gas 
flow  and  propellant.  Therefore,  if  discrepancies  between  measurements 
and  calculations  exist  and  arc  l.argcr  than  can  be  caused  by  the  first 
two  error  sources,  then  one  can  conclude  that  the  mathematical  model 
01  tin;  flo’w  is  inadequate.  In  an  analysis  of  the  response  of  tlic  appa- 
ratus to  theoretical  flows,  errors  from  the  third  source  ui'e  not  exist- 
ent, because  fo’i  theoretical  flows  the  equation  of  state  is  given  by 
definition . 

In  Section  6 of  the  report  we  summarize  the  conclusions  which  can 
be  drawn  from  the  simulated  response  of  the  apparatus  and  give  quantita- 
tive limits  for  application  of  the  measurement  technique. 
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2 . THi;  ACOUSTIC  TIlMPIiRATURn  M»:ASURl:M;iN'T  Tt3aiNIQUl:. 


Technical  details  of  the  acoustic  temperature  measurement  technique 
and  some  tspical  results  are  given  in  Reference  1.  In  this  report  we 
shall  present  only  a short  outline  of  the  method. 

The  apparatus  is  designed  to  measure  uiio  sound  speed  c[m/s]  in  the 
ga.s  behind  the  projectile.  Assuming  that  the  gas  is  adequately  described 
by  an  ideal  gas  equation  one  can  compute  the  gas  temperature  T[K]  from 
the  sound  speed  by 


T = c'M/(YRj 


2 

^ Aflame 


/(sf). 


(1) 


whore  M[kg/mol]  is  the  molar  mass  of  the  gas,  Y is  the  ratio  of  its  spe- 
cific heats,  R = 8.3143IJ  K'^  mol"^]  is  the  universal  gas  constant, 
is  the  isochoric  flame  temperature  of  the  propellant  which 

drives  the  projectile,  g = 9.R0G55[ra/s‘’]  is  the  standard  acceleration, 
and  f[in]  is  the  explosive  "force"  of  the  propellant. 

Tl'.c  sound  speed  c is  obtained  by  subtracting  from  the  projectile’s 
muztlc  velocity  the  velocity  u^,  of  tlie  rarefaction  wave,  w’ljch 

travels  upstream  into  the  barrel  after  the  projectile  has  cleared  the 
mutilo,  1 .c. , by 


(2) 


The  setup  of  the  aj'paratus  for  the  measurement  of  and  is  .shown 

schematically  in  Tigure  1.  The  muztlc  velocity  u^^  is  obtained  from 

X-rays  of  the  projectile  outside  the  weapon.  Tlie  velocity  of  the  rme- 

facrioil  wave  is  determined  from  pressure  traces  wliich  arc  produced  by  two 
or  more  pressure  gages  in  the  barrel. 

l.ct  Ax  = Xn  - X.  be  the  distance  between  the  first  and  the  last 
I)  A 

pressure  gage  (Figure  1).  Then  the  sound  speed  e can  be  computed  by 


c 


(5) 
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if  one  assuinor.  that  the  chai-acteri stic  AUM  in  Tijiurc  1 is  straight  and 
the  particle  velocity  is  equal  to  tliroughont  the  iwi7,£lo  region. 

Typical  pressure  histories  at  the  gages  ai'e  shown  in  I'ig.  2,  They 
were  computed  using  a mathematical  model  [2,3]  of  the  interior  ballis- 
tics and  they  agree  qualitatively  with  pressure  traces  actually  observed 
and  reproduced  in  Reference  1.  The  arrival  of  the  rarefaction  wave  cor- 
ro.s]ioiids  to  a gradient  discontinuity  in  tlic  pressure  trace.s  and  the  arri- 
val times  can  be  determined  by  visual  inspection  of  photographs  of  the 
traces . 

The  simple  formula  (3)  for  the  sound  speed  is  exact  if  all  flow 
properties  within  tlie  area  ABMDC  of  I'igiirc  3 are  uniform.  Making  oflier 
assumptions  about  the  flow,  one  can  derive  ottier  formulas  for  the  sound 
speed.  Tims,  assuming  in  ABMl'C  isentropy  and  adiabatic  expansion,  but 
permitting  the  press\n-c  to  be  variable,  the  following  corrected  formula 
is  derived  in  Reference  1; 


c.,  = 


('  ■ ’N 


r 1 

1 1 

Y - 1 

J 

>-  t . 

-1 


In  eq.  (4),  and  Pj^  arc  the  sound  speed  and  pressure,  respectively, 

at  node  i'  oi  riguic-  3,  and  .'it  = - tj,.  The  integral  in  oq.  (4)  is 

a lino  integral  along  the  path  llA  of  the  raref.ict ion  v\'avc.  11  the  pres- 
sure p 13  pp  along  the  path,  then  eq.  (4)  reduces  to  eq.  (3).  One  can 

sec  from  rigurc  2 that  generally  p/pp  1 along  the  integration  path 

and,  tlmroforo,  c^  computed  with  oq.  (4)  is  smaller  tlian  tl\c  sound  speed 

computed  with  oq.  (.3).  The  integral  in  eq.  (4)  can  he  evaluated  numeri- 
cally using  pressures  which  are  provided  by  pressure  gages  between  tbe 
extreme  gages.  The  appat'Otus  described  in  Reference  1 had  throe  gages, 
located  at  25mm,  fiOmm  and  75min  from  the  muccle  of  a 20mm  cannon.  This 
luimhcr  of  gages  is  sufficient  for  the  numerical  quadrature  because  the 
pressure  varies  smoothly  aiuug  t'lc  integration  patli  and,  in.  fiddition, 
the  pressure  exponent  (y  - l)/dY  is  very  small  (ty])ically  about  0.1). 
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"A,  K.  R.  CclminS,  "TlKmrctical  Basis  of  the  Recoillei.s  Rifle  Interior 
Ballistics  Code  RhCRIF",  BRL  Report  No.  1931,  September  1976.  CAD  ..H(il 38.32!,} 

A,  K.  R,  Cclmip.^,  "RHCHIT  Users'  Manual",  BRl-  Memorandum  Report  2692, 
October  1970. 
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3.  STANDARD  OBST.RVATIONAI,  HRRORS. 


The  principal  idea  oT  the  acoustic  temperature  measurement  techiiiquc 
is  to  measure  the  sound  speed  c and  to  compute  the  corresponding  tem- 
perature T by  cq.  (1).  Because  T is  proportional  to  the  relative 
standard  error  of  T is  twice  the  relative  standard  error  of  c; 

e^/T  = 2 , (5) 

{■q,  (S)  provides  an  estimate  of  from  estimates  of  c^. 

The  standard  error  of  the  sound  speed  depends  on  the  obser- 
vational errors  as  well  as  on  the  formula  wliich  is  used  to  compute  c. 
The  simple  formula  (T)  has  a smaller  standard  error  than  eq.  (4)  because, 
it  depends  on  lesser  measurements.  The  difference  bctn'cen  the  accuracies 
of  botli  equations  is  small,  liowever,  because  the  additional  observations 
(pressure)  enter  eq.  (4)  witli  a very  small  exponent.  We  shall,  there- 
fore, restrict  our  analysis  to  eq.  (3). 

Tor  the  purpose  of  error  analysis  one  can  consider  the  following 
quaiuitics  as  directly  and  independently  observable:  the  muzzle  veloc- 
ity the  arrival  times  of  the  rarefaction  wave,  t^  and  tg,  and  the 

distance  between  the  pressure  gages.  The  law  of  error  propagation, 
(Reference  4)  applied  to  cq.  (.3),  then  furnishes  a formula  for  the  stan- 
dard error  of  the  sound  speed  in  terms  of  the  standard  errors  e ,,,  e 
and  0 of  the  observables.  The  formula  is  ^ 

AX 


( /.  \ 2r/c.  \2  /e.\2-i 

. = <(~1  + 2 i + 

c I ^ / |_\  \ 


wIktc  At  = - tj^!  The  relative  standard  error  of  ‘Jte  sound  speed  is 


2 . .2 


■>  7 


In  order  to  get  an  idea  about  the  relative  importance  of  the  various 
terms  in  cq,  (7),  we  substitute  t)qiical  values  into  the  equation.  The 
ratio  is  tlic  Mach  number  of  the  flow.  In  the  examples  in  this 

report  it  varies  between  0.2  and  0.8. 


W,  I-..  Demming,  ^t^tj_stical_  Adjustment  of  Data,  J.  Wiley  8 Sons,  New 
York,  NY,  1944. 


j 


111  practical  applications  the  Mach  number  will  be  limited  probably  to 
values  below  0.6  because  for  larger  values  the  rarefaction  wave  is  dif- 
ficult to  detect  from  the  pressure  traces.  As  a convcnioiu  ty^iica]  ave- 
rage one  might  assume  Uj^/c  *=  0.5.  The  values  of  and  we 

assume  to  be  of  the  order  0.01.  With  these  assumptio  s eq.  (7)  simpli- 
fies to  ' 


The  standard  error  of  the  time  readings  increases  with  the  muz- 
zle velocity  because  the  rarefaction  wave  is  less  pronounced  in  pressure 
tracc.s  for  higher  velocities.  The  deterioration  of  time  readings  is  com- 
pensated to  some  extent  by  an  increase  of  At  with  increasing  u.,.  For  a 

Mach  number  w 0.3  one  finds  from  Figure  7 of  Rcfercnco  1 that  At  s;  0.1 
ms  and  c^/At  0.1  or  larger.  This  shows  that  in  oq.  (7)  or  (8)  the 

y-  time  reading  inaccuracies  dominate  over  all  other  ,erms.  The  relative 
standard  error  of  the  sound  speed  is  with  these  assumptions  about  7%. 
The  corresponding  relative  standard  error  of  the  temperature  is  about 
14 0.  Thi.s  value  agrees  in  the  order  of  magnitude  with  the  scatter  of 
about  IQo  of  observed  temperatures  shown  in  Figure  12  of  Reference  1, 
Our  analysis  shows  that  a scatter  of  this  magnitude  can  be  ex- 
plained by  inaccuracies  uf  time  readings.  It  is  also  evi''nt  that,  in 
order  to  improve  the  overall  accuracy  of  the  tu<  imique,  one  has  to  im 
prove  in  the  fir.st  place  the  accui'acies  of  the  rarefaction  wave's  arri- 
val time  measurements. 

The  above  consideration  is  restricted  to  random  observational  errors. 
In  addition  to  these,  one  ha.s  systematic  errors,  intrinsic  to  the  tech- 
niquo.  We  shall  analyze  some  of  the  systematic  eners  in  the  next 
sect  ions . 


4.  TllbORI-.TICAL  TIiST  CASllS. 

Reference  1 describes  tests  of  the  acoustic  temperature  mcasuremonL 
technique  on  a 20mm  cannon.  The  tests  were  carried  out  lor  muzzle  veloc- 
ities between  265  m/s  and  910  m/s.  The  variation  of  the  muzzle  velocity 
was  achieved  by  a variation  of  the  amount  of  propellant  which  was  used 
for  the  firing  of  the  weapon. 

In  this  report  wc  chose  an  identical  arrangement  for  the  generation 
of  theoretical  test  cases.  The  interior  ballistics  were  computed  using 
tlic  computer  code  RF.CRIF  [2, .3].  The  code  was  slightly  modified  for  tlic 
present  purposes  to  allow  an  explicit  tracing  of  tlic  rarefaction  wave 
and  to  permit  a better  definition  of  the  flow  within  the  expansion  fan 
region.  The  input  parameters  for  the  calculations  arc  listed  in  Table  I, 
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^ ^ z iin  == ^ ^ 
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Projectile 

>la ;;  s 
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:--rr-r:jj^Q{,  KK/'"’ 
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2(.oo  i; 
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propel  lain 
■■“Density 

Halio  ol  siu'i’iCio  iu-ats  y 
J liinic  lciii]if!'8t,ut'e 
I,x])k':ilvv  "J'oi'cc" 

Cyiiitniiili'jii  lnw 
Cciiiil)Ui;li on  lav;  covlliclvnt 
(icoinvliy:  sjiliv'res  inl'i  Jiunieli;! 

l.tj'iaiion  ol  slate;  Ideal  {jaa  wltli  coiisliinl  if 
Mays:  3i',  - 47, M;  (sevcial  cases  well’  cwij'iiU'd; 

I’lessiiie 

Si  X iMi'i'.s  ill  ir.rniin  jnlui  vals  from  llie  miuxlo.  I oi  slinnliUcd 
tcn.HMalu.e  ineasurenieitls  n ’'rii-.t  KiiKo"  will,  a:, Mimed  ul  2:.mm  un-l 
..  -I  ,.i  ,i„.  iii.-aiion.  Joi  ud)ii;ilmuiil  colcuiul lyii*- 

41  

ull  six  j'.ajp'S  were  used, 
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Iiy  Hjifcl  1 /ill};  clHli-runl  amounl!;  of  j>roiicl]aiit.,mu2z;li;  velocities  between 
2l)U  ni/s  (iiiill2()()  ni/s  wele  obl&lneJ.  The  relation  Itetween  jtioiie)  1 ant  mass 
ami  iiiuzrle  velocity  }s  uiiowii  in  figure  3 which  also  shows  exiierimental  ly 

^hservoil  velocities,  because  tiic  purpose  of  llie  eoltulal  ions  was  to  gen- 

' ^ "crate  rcpTcaentntive  theoretical  interior  ballistic  flows  (and  not  to 

model  a S|iecillt'  wcajioii),  the  agreenent  between  experimental  and  calcu- 
lated velocities  is  sufricicntly  close.  Since  calculated  and  observed 
pres'iUies  also  agieed  uppi oximately , the  calculated  flows  cun  be  con- 
sldeied  (is  uilequaie  exuiiijiles  of  tyjileal  tlicorel j cal ly  computed  flows  In 
(I  20!ani  cuiiiioii , 

In  the  next  section  we  nhall  describe  the  si  mu  la  ted  resi'Onse  of  the 
acoustic  Icinpenitiii  e measurement  ai'paialus  to  the  tlicorei  leal  ly  coinjiuted 
ll<<ws.  To  ttlimilalf  this  response,  one  has  to  extract  from  the 
tulculai  ions  tin  ni  rival  times  Of  the  latefaclion  wave  .it  };a);c  locations, 
list  the  piesSoie  lilstorles  Ot  the  same  locutioiiH,  and  use  these  data 
ill  the  t tiiii’ei  at  111  e f oi  mill  as.  because  tlio  intorloi  ballistics  calcula- 
tions Included  an  explicit  tracing  of  tlie  J ai ef act  ion  wave  front, 
the  anival  tunes  wei'e  obtained  wltli  higb  accuracy,  This  cl  1ml- 

Jiatuil  the  m.ijoi  Koui  ce  of  obset  vat  lonnl  erjot  s.  A secund  source  of 

■”  ■ ■■■ -error,  the  Innccutacy  of  the  eqimtlon  of  *tote,  is  not  prc-sent  in  a »lm- 

lilatc'd  tcsjionse  to  ii  Iheoreticul  flow,  because  the  equation  of  atatc 
is  given  1/y  <li  f lull  Ion,  fgusequei'lly,  any  di  1 fei  ences  lielwcen  the  };as 
t eiiipf'i atui e*  fioiii  Intel  lor  builistlcH  cnlculm  Ions  and  from  the  uiiiiul't 
ltd  response  ol  the  acoustic  nieajUtemonl  technique  aie  caused  by  avei - 
(iglii);  elteci:i  ol  I lie  technique. 

Ibe  av(ia};lii}i  effects  me  umall  if  llit  gas  lempeiatuie  vuiies  only 
: lllll'  wltblii  the  region  AliMbC  of  I'lgute  ],  'Ihe  vai  iaiions  cdiuioi  be 

, (issumed  Sliiiill  in  all  cases,  huweVei  , A one  dimtnsionn I lieutiiient  of  In- 

I Icil'ji  I'.'illlstlcN  produce!,  iieccssai  1 ly  a lomperatui e peal  at  the  basi/i 

of  llie  Jit sjJ ect  1 1 e I'l,*').  Till'  iheoicllcol  baci.gi ouiid  of  the  pheiiomenoii 
i is  discussed  In  l!clcreiice  2,  iq'.  lo)  lu.',  aii'l  a typical  example  (,if  ibeo- 

r I'tlially  predicted  leiiipej  ului  c vm  i at  J oris  wiiblii  tlic  guu  is  ihown  in 

! Il,;uic  4,  111''  teiiipei  at  Ui  e variations  wUhlli  the  l'-gioll  Al'''lli!;  of  ol'Sci 

r Valloiis  can  be  quite  lai};"  foi  lilgli  iiiu.;j;le  ve  loci  I i e:i , llils  is  lllus- 

I Hated  by  I Iguie*  S and  (>,  showing  theoi  et  Icii ) temiieraturc  histories  at 

! the  locutions  ol  the  exiicinc  gitgen  in  .s  low  velocity  and  u model  ate 

; vclotH/  cas",  ii'jijicciiveiy.  In  ihv  low  Velocity  case,  llg'ue  5,  tb'c 


‘ft,  SIk'IUiIi,  a.  IWjigli'i  1111(1  p.  H'ihu,  "bludy  ol  Ib.iit  limisfoi  and  liosloii 
bate  ill  (am  lloiiel;i",  All  lorce  Armameiii  boboi  a lor  y 'Icclili  1 ca  1 I'cpoi  t 'll'- 
7,1  O'J,  Mm  i ll  l'.)7.1, 

^b,  IlclsiT,  "llod(j};r ii|ili'.iivtu  1 alii en  fUr  liisiulioiiilje  flndlmcii'iloiialc 
llassniJiiiuiig'Mi",  bill)  'Ihesls,  'lochiilcul  Ur.veisily  Auclicn,  llb.l,  lieccmbei 
lb'/:., 
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temperature  difference  between  the  nodes  C and  A is  60K  or  5%,  i.c., 
well  below  the  standard  observational  error.  In  the  case  shown  in  Fig- 
ure 6 the  difference  is  about  lOOOK  or  30%  of  the  initial  muzzle  flow 
- -:^iic:ctemperature . One  can  expect  a strong  averaging  effect  by  the  acoustic 
measurement  technique  in  this  case.  r-  --  j. 

The  temperature  measurement  technique  breaks  down  at  about  900  m/s 
for  the  weapon  type  considered  here.  This  limit  is  reported  in  Ref- 
erence 1 and  it  is  caused  by  difficulties  to  detect  the  rarefaction  wave 
in  pressure  traces  from  high  velocity  firings.  Figure  7 shows  a theo- 
retical example  of  such  a pressure  trace.  The  gradient  discontinuity 
at  the  node  B in  Figure  7 is  barely  visible.  The  reduction  of  the  dis- 
continuity is  caused  by  changes  of  the  characteristics'  slope  with 
increasing  flow  velocity.  Figure  8 shows  characteristic  curves  in  a 
case  with  moderate  muzzle  velocity  and  Figure  9 shows  corresponding 
curves  for  a high  muzzle  velocity.  .As  the  Mach  number  of  the  flow 
approaches  one,  the  angles  of  intersection  between  the  rarefaction  wave 
(a  (u-c)  characteristic) and  the  gage  lines  x = const,  become  more  acute. 
Consequently,  the  corresponding  gradient  discontinuities  in  the  pressure 
.r-T . become  less  pronounced.  A theoretical  limit  of  the  method  is 

given  by  a still  higher  muzzle  velocity  at  which  the  rarefaction  wave  is 
swept  out  through  the  muzzle  before  it  reaches  the  pressure  gages.  In 
the  present  example  this  occurs  for  a muzzle  velocity  of  about  1160  m/s. 


Wc  reiterate  that  the  described  flow  properties  arc  based  on  one- 
— dimensional  mathematical  modeling  of  the  flow.  Experimental  verifi- 
cations of  the  d;'tails  do  not  exist  and  one  application  of  the  acoustic 
temperature  measurements  could  be  to  provide  such  a verification  or 
a rejection  of  the  theory. 

S.  SIMUIAT!;!)  RESPONSE  OF  THE  APPARATUS  TO  TllEORPTICAL  MUZZEE  FLOWS. 


The  response  of  an  acoustic  temperature  measurement  apparatus  to 
theoretically  computed  muzzle  flows  can  be  simulated  as  follows.  First, 
one  obtains  from  the  compiited  flow  accurate  "readings"  of  the  arrival 
times  of  the  rarefaction  wave  at  the  pressure  gage  locations.  The  time 
readings  are  substituted  in  eq.  (3)  to  compute  the  sound  speed,  and  the 
temperature  is  computed  using  eq.  (1).  The  corrected  sound  speed  which 
is  given  b>  eq.  (4)  can  be  computed  using  pressures  along  the  integra- 
tion path,  In  the  present  examples,  the  integral  in  eq.  (4)  was  calcu- 
lated by  a trapezoidal  formula  using  pressure  readings  at  six  pressure 
gage  locations  at  12.5mm  intervals. 

The  results  of  the  calculations  arc  shown -in  Figure  10.  The  extreme 
curves  in  the  figure  are  those  of  the  temperatures  and  T^  at  the 

nodes  I)  and  A of  Figure  1,  respectively.  The  temperature  at  the  muzzle 
node  M is  practically  equal  to  T^.  The  difference  between  Tj^  and 
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is  quite  remarkable  in  cases  with  high  mur.slc  velocity.  I'or  mur.zlo 
velocities  over  600  m/s,  it  is  larger  than  lOOOK.  The  acoustic  tem- 
perature measurement  technique  produces  a temperature  between  and 

as  expected.  Somewhat  unexpected  is  the  large  effect  of  the 

isontropic  correction  of  the  temperature  by  eq.  (<1)-  In  the  present 
example;-,  the  temperature  computed  with  cq.  (4)  is  much  closer  to  than 

to  The  intent  of  the  measurement  is  to  ol.itain  a value  for 

(or  . Obviously,  the  .simpler  eq.  (.i)  is  preferable  to  cq.  (4)  for 

'this  purpose.  The  reason  for  the  failure  of  cq.  (4)  to  improve  the 
tempo’-aturo  is  that  it  is  based  on  assumptions  about  flow  properties 
winch  are  not  met  by  the  computed  flow,  i.c.,  by  the  solution  of  the 
one-dimensional  flow  equations. 

figure  10  shows  tliat  the  acoustic  measurement  technique  based  on 
eq.  (3)  is  applicable  to  flows  with  muzzle  velocity  of  about  SOO 

m/s  or  less.  At  SOO  m/s  the  systematic  difference  between  and  the 

measured  temperature  approaches  the  magnitude  of  the  observational 
errors.  An  extension  cf  the  method  to  cases  witli  higher  muzzle  veloci- 
ties requires  a better  definition  of  the  rarefaction  wave  front.  If 
cq.  (3]  is  used,  then  one  makes  the  implicit  a.'sumption  that  the  wave 
front  is  a straight  line  in  the  x,t-plane.  Tic  inverse  slope  of  that 
line  is  used  in  cq,  (.t)  as  an  approdsiacion  of  the  rarefaction  wave 
velocity.  An  inspection  of  Figures  3 and  9 suggests  that  the  wave  front 
Can  be  better  approxiuUited  by  a third  degree  poUntomial  in  t.  However, 
in  order  to  determine  such  a curve  one  needs  more  than  the  three  nodes 
corresponding  to  the  three  pressure  gages  of  the  experiment.  Thoreforc, 
we  assume  that  tlie  device  consists  of  an  array  of  six  gages,  located  at 
12.5mm  intervals.  (Such  a device  has  been  also  tested  at  URL.)  Then 
tlic  task  is  to  fit  a third  deg^'oe  pol)niomial  to  six  nodes  or  to  seven 
nodes  if  the  muzzle  time  is  included  as  a seventh  observation. 

Let  X.  be  tlic  distances  of  the  six  pressure  ga;,es  from  the  muzzle 
and  t.  bo  ^thc  corresponding  arrival  times  of  the  rarefaction  wave.  Tl.cn 
the  constraints  for  the  curve  fitting  piohltm  are 


(H 


i = 1 , 2 , . . . , 6 . 


The  problem  has 

zlo  time  is  used 
i s s i rap  1 y 


four  parameters,  namely  and  U,^.  If  the  muz- 

as  a seventh  observation  then  '..he  constraint  foi-  it 
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After  the  adjustment  the  rarefaction  wave  front  x (t)  and  its  vcloc 
ity  u^,(t)  can  be  computed  by  ■ ~ — ; --- 

*r<')  = ^l  • “m"^V  * 

and 

“r<'>  ■ (12) 


At  Iho  muzzle  the  wave  has  the  velocity 


- ^M> 


-^and^  t corresponding  temperature  is 


(13) 


T 

■j 

«•- 
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One  can  also  compute  the  temperatuz'o  atong  the  rarefaction  wave  by 


\M  ‘ (uct) 

particle 


(15) 


provided  that  the  particle  velocity  is  known.  That  velocity 

is  generally  not  constant  in  the  muzzle  region,  as  figure  il  s’  owe. 
Only  in  cases  witli  low  muzzle  velocities  the  particle  velocity  can  be 
ap])roxiinated  by  the  muzzle  velocity  Uj^.  These  arc  the  same  cases  where 


llic  lompcraturc  is 
Tlici'efore,  eq.  (15) 


approximately  constant  within  the  muzzle  region, 
is  useless  for  practical  applications  and  wc 


restrict  ourselves  to  the  computation  of  the  muzzle  temperature  T,  by 
cq.  (14J.  ” 


The  data  fittings  were  carried  out  using  the  theoretical  data  and 
assuroii.g  the  same  estimates  of  observational  errors  as  in  Section  3, 
name  1 y 


Cj,  = 0.01  ms, 
c =0.1  mni, 

X 


(16) 


27 


t 


liach  fitting  produced  a value  of  v/liich  then  wa:.  uricd  in  cq.  (Hj  to 

compute  The  results  are  sliown  in  figure  10,  'Mu-  approximation  is 

.^excellent  up  to  muztle  velocities  of  900  in/s,  l.c,,  up  to  the  limit 
where  rarefaction  wove  arrivals  con  be  practically  oluiervcd. 

In  spite  of  the  good  results  in  the  llicoretical  lest,  a lliijd  degree 
iip])roxim;U  ion  of  the  wave  frofit  should  he  used  Judiciously,  liccaose  0) 
the  sensi t i vi ty  of  (he  results  to  ohserval iona i errors.  The  alandurd 
errors  of  the  fittini;  parameter  U,,  and  oi  1 can  be  coiiijnited  using 

variance  propagation  foitnulas  and  assuming  data  st.indard  errois  as  given 
by  eq.  (lo).  The  results  aie  shov/n  in  liguie  IJ,  Tiiey  suggest  that  for 
low  inuctle  velocities  a tbiid  degree  ai)proxi  m.'ii  i on  induces  an  unaccej)!- 
ahly  large  uncertainty  of  the  le(n])et aiure . 'Ilu'  reason  fot  the  lirgi 
vaiucs  of  estimated  standard  ei  rors  is  the  sinall  cuivature  ul  the  jnre- 
faction  wave  liont,  hy  fitting  a higli  order  polynoinlal  to  .sltalghi  line 
data  one  introduces  redundant  paraiiieicrs , but  paramet  ei ::  o)  leduudant 
sets  a-re  strongly  conclated  and  tyi)icalJy  bave  large  Indjvidual  c)  tor 
estimales.  'llie  error  estimates  of  ftom  a lineai  fit  (wluctt  doett  not 

have  redundant  parameters)  nrc  mucli  smaller  and  of  the  same  order  as 
derived  in  Section  3.  Theiefore,  one  should  test) let  ihlid  ordet  fil- 
lings ol'  tlie  rirelaclion  wave  front  lo  cases  whole  the  lionl  i;>  clcaily 
tiOl  a st/aighi  line, 

'file  inclusion  of  the  iiiu;'.r.le  lime  ar.  a seveiitli  oliservali on  in  the 
adjustineiit  problem  docs  not  change  the  lomjicrature  sign!  I ioaiilJy . ilu' 
seventli  olistrvaiion  docs,  liuwcvcj , i educe  the  esliiimtes  ol  llie  standard 
enors, 

C,  CON'tibllSKJNS. 

llie  acoustic  icmpct aiure  measurement  ajipaialir'  niea  uies  an  av'eiagir 
lemjierature  tbi  the  tiisl  I ew  mi  11  i seconds  in  the  mu 1 e flow,  'Ihe 
standard  error  ol  th''  measured  temj"Talure  Js  about  ll>», 

One  has  the  lullowing  guideline:,  for  llie  compiiiisoii  of  expei  Itiieni  ai  ly 
measured  lemiierat ui es  with  tlieorelieal  taieuliilions  ol  the  same  events: 

o Hi  sagreement  within  1!.>1  is  loJeiabJe, 

V 1 eiiipei  atU)  e peaks  ol  less  ihan  one  ml ) 1 ; :ieconil  duiatjoii,  II  exist- 
cm  , cannot  be  delected  by  the  aj>pai;itus,  because  ol  its  uveiaging 
property. 

0 lii  sugreeiiienls  exceeding  leT  In  a I low  wllliout  shoi  t duialJon  lem- 
jiciatuie  peaks  indieute  inadequate  inydelliig  ol  the  flow, 

li  llie  acoustlculjy  meo'iuied  t eiri]>t)i  atu)  o Is  used  Lo  i baj  iu.  Ici  1 ;.e  tlw 
f low,  CM  Js  coinji.’ii ed  with  tcmpoiature  moaswi omeiit :<  1/  olhel  loebnJ'iuea, 
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A.-u.acy  or  Ui.laiKfJ  !->'  Unvt 

liUiiiK  I'uiclacl  iuit  Vtiivc. 


,Ke„  one  should  keep  in  .ind  the  po.slhility  o£  a syato..tic  error  due 
to  an  approximate  equation  of  state. 

The  standard  error  »£  the  .eaeupd  Sf  r”!'. 

niflcantly  by  increaaing  the  accur.ei-  ^ reduction  of  standard  observa- 
Jronlrcro'rs'  r/'So  'te:;:;ature  could  be  achieved  by  ..aKing  .easure- 
inents  of  repeated  identical  firings. 

• cv,  T effect  of  the  technique  is  possible  by 

A reduction  ot  vue  avciugii.*  ..c-im,  a hicher  order  apt<roxi- 

increasing  the  number  of  pressure  8^2°-  approach  should  be  used 

Nation  for  the  'Te';itIJ  y of  the  te.petaturo  to 

judiciously  because  it  incicases  the  sensitivny  r 

ol)scrvational  errors. 
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